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ABSTRACT 

We present observations of the unusually bright and long 7-ray burst GRB 050820A, one of the best-sampled 
broadband data sets in the Swift era. The 7-ray light curve is marked by a soft precursor pulse some 200 s before 
the main event; the lack of any intervening emission suggests that it is due to a physical mechanism distinct 
from the GRB itself. The large time lag between the precursor and the main emission enabled simultaneous 
observations in the 7-ray, X-ray, and optical band-passes, something only achieved for a handful of events to 
date. While the contemporaneous X-rays are the low-energy tail of the prompt emission, the optical does not 
directly track the 7-ray flux. Instead, the early-time optical data appear mostly consistent with the forward 
shock synchrotron peak passing through the optical, and are therefore likely the beginning of the afterglow. On 
hour time scales after the burst, the X-ray and optical light curves are inconsistent with an adiabatic expansion of 
the shock into the surrounding region, but rather indicate that there is a period of energy injection. Observations 
at late times allow us to constrain the collimation angle of the relativistic outflow to 6.8° < 9 < 9.3°. Our 
estimates of both the kinetic energy of the afterglow (£ke = 5-2^-j x 10 51 ergs) and the prompt 7-ray energy 
release (E 1 = 7.5^ x 10 51 ergs) make GRB 050820A one of the most energetic events for which such values 
could be determined. 

Subject headings: gamma rays: bursts — X-rays: individual (GRB 050820A) 



1. INTRODUCTION 

With the discovery of the cosmological n ature of 7-ray 
bursts (GRBs) in 1997 JMetzger et all 119971) . astronomers 
were suddenly forced to explain the enormous isotropic en- 
ergy release of these distant explosions. Some of the most 
energetic events, such as GRB 990123, seemingly released 
enough en ergy in the prompt 7-rays (£^0 = 1.2 x 10 54 ergs; 
Brigg s et alJl99 9) to rival the rest mass of a neutron star. Fur- 
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thermore, broadband modeling of the best sampled events has 
shown that a comparable amount of energy remains in the 
shock, powerin g the long-lived X-ray, optical, and radio af - 
terglow (see e.g. Panaitescu & Kumar 2001; Yost et al. 2003). 

The hypoth esis that GRBs are a-spherical explosions 
jRhoadslll9 99). supported by the appearance of achromatic 
"jet" breaks in a large number of afterglow light curves 
ISari et alJ 1 1999ft . proved to be a turning point. With typi- 
cal opening angles of a few degrees, the true energy release 
from most GRBs is ~ 10 51 ergs, on par with that of a super- 
nova (SN). This realization enabled the discovery of a stan- 
dard en ergy reservoir fo r the collimation-corrected prompt 
energy (|F rail et alJl2001) and kinetic energy of the afterglow 
( Berger et al. 2003a). GRBs are now considered promising 
standard candle candidates, with the hope of Hubble diagrams 
out to z ~ 6 offering complementary c onstraints to Type la 
SNe on the cosmology of our uni verse (Firmaniet al. 2006; 
iDai et alJ2004l c.f. lFriedman & BloorrJ2005l) . 

Launche d in Novembe r 2004, the Swift Gamma-Ray Burst 
Explorer ( Gehrels et al. 2004J) was designed to position 
GRBs, disseminate accurate coordinates to ground-based ob- 
servatories in real-time, and follow the UV and X-ray after- 
glows from minutes to days after the event. In only a year 
of full operation, Swift has brought about a number of fun- 
damental advances in the GRB field, inclu ding the d iscov- 
ery of the first X-ray (GRB 050509b:lGehrels et all200l and 
near-infrared (GRB 050724: lBerger et all2005l) afterglows of 
a short-hard burst, the detection of the high-redsh ift (z = 6.3) 
GRB 050904 (Haislip et a lJ2006tlCusumano et alJ2006h . and 
the ability to measure broadband light curves starting shortly 
after, and in a few cases even during the 7-ray event itself. 

Despite these advances, measuring the bolometric fluences 
of Swift events has proved challenging, for a number of rea- 
sons. First, the limited energy range of the Swift Burst Alert 



2 



Cenko et al. 



Telescope (BAT; Barthel mv et alJ l2005) means that Swift can 
accurately characterize only the softest GRB spectra. Sec- 
ond, few Swift events have shown conclusive signs of a jet 
break, leaving geometric corrections highly uncertain. Fi- 
nally, the X-ray ligh t curves of Swift aft erglows have shown 
both bright flares llBurrows et al.l l2005bl) and slow decays 
JNouseketalJl2 006). Both behaviors have been attributed to 
late-time (t >> ?grb) energy injection, and at times have ri- 
val ed the energy releas e of the prompt 7-ray emission (see 
e.g. Falcone et al. 2006). 

On 20 August 2005 UT, the BAT detected and localized the 
unusually bright and long GRB 050820A, a truly rare burst in 
the Swift sample. The 7-ray light curve is marked by a soft 
pulse of emission preceding the main event by over 200 s. 
The main emission was bright enough to be detected by the 
Konus-Wind instrument, providing a 7-ray spectrum extend- 
ing beyond 1 MeV as well as continuous coverage over the 
entire ~ 600 s burst duration. 

Since Swift triggered on the precursor, both space- and 
ground-based facilities were able to image the transient dur- 
ing the bulk of the prompt emission. Such contemporaneous 
multi-wavelength observations have only been achieved for a 
handful of bursts to date. The bright X-ray (F v ~ 0.7 mJy) and 
optical (R ~ 14.5 mag) afterglows made it possible to study 
the evolution of the afterglow for weeks after the burst, pro- 
viding one of the most detailed broadband light curves in the 
Swift era. Finally, late-time Hubble Space Telescope obser- 
vations allowed us to constrain the jet-break time, and hence 
the geometry of the outflow. Even after applying the colli- 
mation correction, we find GRB 050820A is an exceptionally 
energetic event. 

Our work proceeds as follows: in 5j2j we outline our broad- 
band observations of the GRB 050820A, beginning with the 
high-energy prompt emission and followed by the X-ray, opti- 
cal, and radio afterglow. We find the afterglow data are incom- 
patible with the standard model of synchrotron radiation from 
a single, highly -relativisti c shock ex panding adiabaticallyinto 
the surrounding medium (ISari et alJll998l) . Instead in |5] we 
use power-law fits (F„ oc ?" C V"' 3 ) to model the afterglow, di- 
viding the burst into segments based on noticeable temporal 
breaks in the X-ray and optical light curves. This anal ysis 
allo ws u s to investigate the early broadband light curve ( q4.ll 
and !4.2>, late -time (t » ?g«b) energy injection in the for ward 
shock fll4.3t . the structure of the circum-burst medi um ( 84.41 . 
and the geometry and energetics of the event ( SI4.5> . 

Throughout this work we adopt a standard cosmology with 
H = 71 km s" 1 Mpc" 1 , fi M = 0.73, and ft A = 0.27. We also 
make use of the notation Qx = 10* x Q. 

2. OBSERVATIONS AND DATA REDUCTION 

In this section, we present our broadband observations of 
GRB 050820A, which span the spectral range from 7-rays 
to radio frequencies and extend in time out to 61 days after 
the burst. We include an independent analysis of the Swift 
BAT data set, as well as the complete light curve and spec - 
trum from the Konus-Wind instrument ( Aptekar et al. 1995), 
which, unlike the BAT, was able to observe GRB 050820A 
over its entire duration (" 32. 1> . In 32.2lwe provide an analy- 
sis of the Swift X-Ray Telescope (XRT; Burrows et al. 2005a) 
data, with afterglow detections out to two weeks after the 
event. We present contemporaneous optical data from the au- 
tomated Palomar 60-inch telescope (P60; Cenko et al. 2006, 
in prepa ration) and the Swif t Ultra- Violet Optical Telescope 
(UVOT; Roming et al. 2005), supplemented by late-time im- 



TABLE 1 

Spectral Properties of BAT 7-RAY 
Emission 



Time Interval 


r 


X ?/d.o.f. 


('bat, s) 






-17 - 22 


1.74 ±0.08 


1.07/75 


217-241 


1.07 ±0.06 


0.95 / 76 



NOTE. — Spectra were fit to a power-law 
model of the form dN/dE oc E~ r . Errors re- 
ported are 90% confidence limits. 



ages taken with the 9.2-mHobby-Ebe rly Te lescope (HET) and 
the Hubble Space Telescope (HST) ( SI2.3> . Finally, we mon- 
itored GRB 050820A in the radio with the Very Large Array 
beginning only hours afte r the burst and continuing for ap- 
proximately two months ( 82.4> . 

2.1. Gamma-ray Observations 
2.1.1. Swift Burst Alert Telescope 

At 06:34:53 on 20 August 2005 UT 17 , the BA T triggered 
and loc ated GRB 050820A (Swift trigger 151207. |Pageet"aT] 
2005b). The initial location calculated on-board was a 4' 
error circle centered at a = 22 h 29 m 35."9, 6 = +19°ll'14f'2 
(12000.01 ICummings et all d2005a) describe a multi-peaked 
light curve (tgo = 26 ± 2 s) with clear spectral evolution (hard- 
to-soft) within each peak. 

Following the report of additional high- energy emission 
from the Konus-Wind instrument shin & Frederiks l2005l 
see ^2. 1 .21 . the BAT team re-analyzed their full light curve 
and found evidence of a much stron ger, harder episo de of 
emission from GRB 050820A ('Cu mmings et al.ll2005bl) . Un- 
fortunately the satellite entered the South Atlantic Anomaly 
(SAA) approximately 240 seconds after the burst trigger, and 
thus estimates of the properties of this second phase are highly 
uncertain. 

Here we have independently analyzed the BAT data from 
GRB050820A. We have extracted the 15-350 keV light 
curve in 1 s timing bins using software tools from the Swift 
data analysis package 18 . The result is shown in Figure^. 
In addition, we have extracted spectra for the two periods of 
high en ergy e mission covered by the BAT data (peaks A and 
B, see ^2l 1 .21 . We then fit these spectra to a power-law dis- 
tribution of energies (dN /dE oc E~ r ). We find evidence for 
strong spectral evolution between these two intervals, as the 
second peak is significantly harder than the first. The results 
of this analysis are shown in Tabled 

2.1.2. Konus-Wind 

The main part of GRB 050820A triggered Konus-Wind at 
7kw = 06:39: 14.5 12 UT, 257.948 s after the BAT trigger (tak- 
ing into account the 3.564 s propagation delay from Swift to 
Wind). It was detected by the S2 detector which observes the 

17 It is customary to refer to the burst trigger time as 7b, for it is assumed 
to coincide with the beginning of the prompt emission. Given the unique 
nature of the high-energy emission from GRB 050820A, we undertake a more 
det ailed study to determine exactly when the prompt emission began (i.e. 7b) 
in i|4. II Times measured with reference to the Swift trigger time will be 
referred to as /bat throughout the remainder of this work. 

18 Part of NASA's High Energy Astrophysics Software package, see 
http://heasarc.gsfc.nasa.gov/docs/software/lheasoft 
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FIG. 1 . — Early Broadband Emission from GRB 050820A. (a): Swift-BAT light curve extracted from 15—350 keV in 1 s bins. The dashed vertical line is the time 
of Swift trigger, 06:34:53 on 20 August 2005 UT (Tbat)- While the second, brighter period of emission is clearly visible, Swift entered the SAA approximately 
240 s after Tbat, effectively terminating the observations, (b): Konus-Wind light curve extracted from 18- 1 150 keV. The 5 peaks visible in the Konus-Wind light 
curve are labeled A — E, and defined in Tablel2l The left dashed vertical line shows the Swift trigger time, while the right dashed line shows the Konus-Wind trigger 
time, Tkw, 258 s later. The portion of the light curve covered by the BAT comprises only a small fraction of the total 7-ray emission, (c): Contemporaneous 
Swift XRT observations (black diamonds) overlaid on the BAT light curve. The X-ray data nicely track the 7-ray emission, (d): Contemporaneous UVOT and 
P60 optical data overlaid on the Konus-Wind light curve. Unlike the X-ray, the optical is not a good trace of 7-ray emission. 



north ecliptic hemisphere; the incident angle was 63. °2. Count 
rates are continuously recorded by Konus-Wind in three en- 
ergy bands: Gl(18-70 keV), G2(70-300 keV), and G3(300- 
1 150 keV). Data collected in this "waiting mode" are acquired 
in 2.944 s timing bins. The time history recorded in the three 
energy ranges can be considered a continuous three-channel 
spectrum. 

Immediately following the Konus-Wind trigger, the instru- 
ment began simultaneously collecting data in "trigger" mode, 
as well. From 7kw to 7kw + 491.776 s, 64 spectra, each 
composed of 101 energy channels ranging from 18 keV - 14 
MeV, were accumulated. The time resolution of these "trig- 
ger mode" spectra varies from 64 ms - 8.192 s and is deter- 
mined by an automated on-board algorithm based on count 
rate. Data were then processed using standard Konus-Wind 
analysis tools and spectra were fit with XSPEC. 

At least five pulses are evident in the Konus-Wind light 
curve (labeled A- E, see Fig. ]l]p and Table |2j. Peak A, 
which generated the BAT trigger, appears to be a weak pre- 
cursor. Figure |2] shows the three-channel light curve of 
GRB 050820A, as well as the hardness ratios. With the ex- 
ception of the precursor (peak A), the burst shows an overall 
hard-to-soft evolution over the entire burst duration, as well 
as within some of the individual peaks (B and C). 

The spectra of individual pulses are well fit by a cut-off 
power-law model: dN/dE oc E~ a exp(-(2 - a)E / E v ); here a 
is the photon inde x and E p is the peak energy of the v x F v 
spectrum (see e.g. ^4.2> . Fitting the overall Konus-Wind "trig- 
ger" mode spectrum, accumulated from < ?kw < 295 (258 < 
f BA T < 553) 19 in the 18-2000 keV range yields a = 1.41+gff, 
and a peak energy E p = 271^f 9 keV (90% c.L; x r 2 = 0.74 for 
62 d.o.f.). 

However, in examining the full Konus-Wind light curve 

9 After this time only a weak tail is seen in the Gl band up to (bat ~ 730 
s; this tail contains less than 5% of the total burst fluence. 



TABLE 2 

Properties of Konus-Wind 7-Ray Light Curve 



Peak ID 


Time Interval 


Fluence a 


Peak Flux 




(?BAT, s) 


(10" 6 erg cm" 2 ) 


(10~ 7 erg cm -2 s~') 


A 


-4.3-19.3 


2.77 


1.7 


B 


222.4-282.8 


28.7 


13 


C 


397.5-430.2 


10.2 


4.3 


D 


454.8-479.4 


3.20 


1.9 


E 


520.4 - 544.9 


4.93 


2.6 


Total 


-4.3 - 544.9 


c 9 7+li.H 
JZ -'-6.9 





a The fluence was extracted from the 20-1000 keV energy range 
(observer frame) assuming a cut-off power-law spectrum of the form 
&N /dE oc E~ a exp _<2_a ' l£/ ' £ P . Errors reported are 90% confidence lim- 
its (see 312. 1.2l for further details). 



(Fig. \$}p), it is clear the above time interval does not in- 
clude a sizable fraction of the 7-ray emission. To derive 
the spectral parameters of the time-integrated spectrum over 
the main part of the GRB (peaks B—E), we simultaneously 
fit the three-channel Konus-Wind spectrum accumulated from 
-33 < fKW < and the overall multichannel spectrum. We 
find a = 1.12^;1| and a peak energy E p = 367^ keV (90% 
c.L; Xr = 0-99 for 64 d.o.f.). Not surprisingly the peak en- 
ergy increased significantly, as the beginning of Peak B was 
the hardest portion of the entire burst. We consider this 
fit the most accurate estimate of the spectral properties of 
GRB 050820A. 

The fluence and peak flux for each individual episode 
are shown in Table |2 The total fluence received from 
GRB050820A from 20-1000 keV (observer frame) was 
5.27^;^ x 10" 5 ergs cm -2 (90% confidence limit). 

2.2. X-ray Observations 
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FIG. 2. — Spectral evolution of GRB 050820A. The top thi'ee panels show 
the Konus-Wind light curve divided into three energy bands: Gl = 18 — 70 
keV, G2 = 70-300 keV, and G3 = 300- 1050 keV. The bottom two panels 
show the resulting hardness ratios: G2 / Gl and G3 / G2. Background levels 
are indicated with horizontal dashed lines. The vertical dashed lines denote 
the BAT (7bat) and Konus-Wind (7kw) trigger times. With the notable ex- 
ception of the precursor, the burst shows an overall hard-to-soft evolution, 
both over the entire duration and within individual bright peaks (B and C). 



The Swift XRT began observations of GRB 050820A 88 
s after the burst trigger. A bright, fading X-ray source was 
detected by the automated on-board processing routine at 
a = 22 h 29 m 37."8, S = +19°33'32"7 (7" error radius) and re- 
ported in real-time (Pageetal. 2005b). Data was taken in 
Window Timing mode until ~ 1 .29 hr after the burst, when, 
due to the decreased count rate, Photon Counting mode was 
automatically initiated. 

Here we have independently processed and reduced the 
XRT data from GRB 05082 0A (for a previous an alysis of the 
early-time X-ray data, see O'Brien et al . 2006). To reduce 
the X-ray data, we used the software tools available from the 
Swift Data Center. We followed standard reduction steps, ex- 
cept for taking measures to mitigate the effects of pulse pileup 
in our spectral analysis. During the steep rise in the X-ray 
light curve at ?bat ~ 220 s, the XRT countrate exceeds the 
200 counts s" 1 at whi ch pileup can affect the detector response 
jRomano et all2006T) . For this segment, we removed the cen- 
tral two pixels from the point source image. At the begin- 
ning of Photon Counting mode, we removed the central four 
pixels of the photon counting mode image for our spectral 
analysis, again in ord er to mitigate the effect of pulse pileup 
(Vau g haner^il2006l. 

We fit the X-ray spectrum for Phases la, lb, and 2 sepa- 
rately (see ^3] for a full discussion of the division of the X- 
ray light curve into phases). We used an absorbed power- 
law model, and fit both with the column fixed to the Galactic 
value (« H ,Gai = 5.0 x 10 20 cm" 2 ; Dickey & Lockmadll990l 
as well as with the column floating. For Phase la, the best 
fit column is consistent with «H.Gai, and there is little differ- 
ence between fits with the column floating and fixed. For 
both Phases lb and 2 we find acceptable fits with fixed col- 



umn, however the best fit columns are 1.5 ± .23 x 10 21 cm -2 
and 1.3 ± .09 x 10 21 cm" 2 respectively. Therefore we find 
marginal evide nce for ab sorption in excess of the Galactic 
value (c.f. Paae et al. 2005al). 

The results of our analysis of the X-ray light curve of 
GRB050820A are shown in Table For the discussion be- 
low we adopt the spectral fits in Table [3] These were used 
to scale the binned count rates to fluxes and flux density at a 
nominal energy of 5 keV. 

2.3. Optical Observations 
2.3.1. Palomar 60-inch Telescope 

The automated P60 responded to GRB 050820A and be- 
gan a pre-programmed sequence of observations starting 3.43 
minutes after the Swift trigger. The system is equipped with 
an optical CCD with a pixel scale of 0.378" / pixel. Images 
were taken in the Rron R and / and Gunn g and z filters. All 
P60 images are processed with standard IRAF l lTodvll 19861) 
routines by an automated reduction pipeline in real-time. 
Manual inspection of the first images revealed a bright vari- 
able source inside the XRT error circle at a = 22 h 29 m 38."l 1, 
S = +19°33'37."1 (see Figure|5J- This position was pro mptly 
repor ted as the afterglow of GRB 050820A iFox & Cenkol 
120051) . allowing ot hers to obtain high-res olution spectroscopy 
of the afterglow ( Prochaska et al. 2005). We continued to 
monitor the afterglow of GRB 050820A with P60 for the fol- 
lowing 7 nights, until it was too faint for quantitative photom- 
etry. 

Optical photometry of the afterglow was complicated by 
the presence of two nearby objects: one R ~ 20.2 mag star 
located 4.0" south-west of the afterglow, and a fainter R ~ 
21.3 mag object located only 2.9" north-east of the afterglow 
(see Figure [3] right panel). On some nights of poor seeing, 
the full-width at half-maximum of our point-spread function 
(PSF) was larger than 2."0. We have therefore performed 
PSF-matched image subtraction using the CPM technique of 
iGal-Yam et alJ fc004l) on our optical data. Errors were es- 
timated by placing 5 artificial stars with flux equivalent to 
the afterglow in locations with similar background contami- 
nation. In addition, we have also used aperture photometry 
(DAOPHOT) to extract the afterglow flux. On the first night, 
the afterglow was bright enough to be well-detected either in 
single images or short co-additions (< 360 s). For these im- 
ages, both nearby sources were below our detection limit. Re- 
sults from aperture photometry and image subtraction were 
therefore consistent. We quote our aperture photometry re- 
sults for these data, as image defects from imperfect PSF- 
subtraction seemed to artificially inflate these errors. How- 
ever, on subsequent nights the afterglow flux was either near 
or below the level of these nearby objects, and we therefore 
report results from our image subtraction technique. 

Photometric calibratio n was p erformed relative to 20 field 
stars provided in Henden ( 2005 ). Rron R is similar to Cousins 
R (Re) and was treated as identical for photometric calibra- 
tion of these images. Magnitudes from the standard John- 
son/Cousins system we re tra nsformed to Gunn g using the 
empirical relation from lKentl d!985h . We found, however, a 
better correlation between our /-band filter and the Cousins / 
(Ic) filter than from t he transformation to Gunn i provided in 
Thuan & Gunn ( 1976). We therefore use Ic in the remainder 
of this work. For the Gunn z filter , we used the Two Micron 
All Sky Survey (2MASS; Skmtsk ieltd]l20 06) catalog and 
the optical photometry provided bv iHendenl d2005h to inter- 
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TABLE 3 

X-ray Afterglow Spectral and Temporal Fits 



Phase 


.start 

'bat 
(s) 


slop 

'bat 
(s) 


a 


x 2 (a)/d.o.f. 


P 


Xr 2 (/3)/d.o.f. 


la 





217 


2.2 ± 0.3 


1.18/11 


0.90 ±0.09 


0.63 / 349 


lb 


217 


257 






-0.10 ±0.03 


1.02/749 


2 


4.8 x 10 3 


8.7 X 10 4 


0.93 ±0.03 


1.14/29 


1.20 ±0.04 


1.00/770 


3 


8.7 x 10 4 


1.7 x 10 6 


1.25 ±0.07 









NOTE. — We separately fit the X-ray light curves and spectra to power-law models of the form 
F v oc V a and F v oc , respectively. The temporal decay in Phases 2 and 3 were fit jointly as a 
broken power-law, with the break time as a free parameter. In Phase 3, we could not meaningfully 
constrain the spectral index due to the low count rate. All errors quoted are 90% confidence limits. 



polate the spectral energy distributions of 10 sources to the 
Sloan z' bandpass. Typical RMS variations in the calibration 
sources were 0.03 mag in Rq and Iq, 0.04 mag in g, and 0.15 
in z'. 

The results of our P60 observations are shown in Table [8] 
While the expected galactic extinctio n in this d irection is 
small (E(B-V) = 0.044 mag: ISchlegel et al.lll998l) . we have 
incorporated it into our results because of the large wave- 
length range spanned by our observations (E(U—z') = 0.17 
mag). Errors quoted are ler photometric and instrumental er- 
rors summed in quadrature. For this and all other optical data 
in this work, magnitudes a re converted to flux de nsities using 
the zero-points reported in Fukugit aet alJ (Tl995). 

2.3.2. Swift Ultra-violet Optical Telescope 

The Swift UVOT automatically slewed to the BAT location 
and began observations only 80 s after the trigger. However, 
the UVOT also becomes inoperable in the SAA and therefore 
does not cover the period of the main 7-ray emission. 

The Swift team reduced the early-time UVOT data and 
report ed detections in the V, B, U, and UVWl bands 
(Ches ter et al 1 12005). Here we have independently reduced 
the U-, B- , and V- band UVOT data following the recipe 
outlined in iLi et alJ (2006) (see §3.6). As a check, we re- 
calculated the B an d V z ero-points with respect to the field 
stars from Henden (2005[). While our zero -points are consis- 
tent with the ones quoted in lLi et af] (1200 6). we found a much 
larger scatter (~ 0.10 mag vs. 0.01 mag) that could not be 
attributable solely to spread in the field stars. We have there- 
fore incorporated the resulting zero-point errors for these data 
points (as well as a similar value for [/-band). The results of 
these observations are shown in Table|S] 

2.3.3. Hobby-Eberly Telescope 

We triggered target-of-opportunity observations on the 9.2- 
m Hobby-Eberly Telescope beginning on the night of 26 Au- 
gust UT. Observations were taken in both Rq and Ic filters. 
A second reference epoch was taken on 29 August UT. Im- 
age subtraction was performed on the two epochs to remov e 
contamination from nearby sources (as described in 32.3. 1> . 
Photometric calibrat ion was perfo rmed relative to 10 refer- 
ence objects from Henden (2005), and the absolute calibra- 
tion was of similar quality to the P60 dataset. Our results are 
reported in Table|S] 

2.3.4. Hubble Space Telescope 

To better constrain the jet-break time, as well as to in- 
vestigate the properties of the host galaxy of GRB 050820A, 



we triggered our Cycle 14 Hubble Space Telescope program 
(GO-10551; PI: Kulkarni). Using the Wide-Field Camera 
(WFC) channel of the Advanced Camera for Surveys (ACS), 
we imaged the field of GRB 050820A on 29 September 2005 
UT (f BAT ~ 37 d) in the F625W (r'), F115W (/'), and F850LP 
(z') filters (see Figure [3}. A second epoch to study the host 
galaxy is scheduled for July 2006, and these results will be 
reported in a future work. 

The HST data were processed using the multi drizzle 
routine (Fruchter & Hook 2002) in the stsdas IRAF pack- 
age. We used pixf rac= 0.8 and pixscale= 1.0 for the 
drizzling procedure, resulting in a pixel scale of 0.05" / pixel. 
The astrometry on these images was then tied to a P60 co- 
addition of all the /?c-band data taken on 20 August 2005 
(which is itself tied to the USNO-B1 astrometric grid). 

The afterglow is well separated from any nearby objects 
in the field, and so we have followed the recipe for aperture 
photometry from Siriann Tet alJ ( 120051) . As a note of caution, 
however, flux from an underlying host galaxy could affect the 
results reported here. Final values for the late time afterglow 
flux will require image subtraction of any host galaxy contri- 
bution, expected following our July observations. F625W and 
F715W magnitudes were converted to the Rc and Iq band- 
passes using synthetic spectra from Table 22 of Sirianni et al. 
( 2005). The results of our measurements are shown in Table[8] 

2.4. Radio Observations 

In Table |4] we summarize our radio observations of 
GRB 050820A, spanning 0. 1-61 days after the explosion. All 
observations were conducted with the Very Large Array 20 in 
standard continuum mode with a bandwidth of 2 x 50 MHz 
centered at 4.86, 8.46 or 22.5 GHz. The array was in the C 
configuration, with an angular resolution of 2.3". We used 
3C48 (J0137+331) for flux calibration, while J2212+239 and 
J2225+213 were used to monitor phase. Data were reduced 
using standard packages within the Astronomical Image Pro- 
cessing System (AIPS). 

3. ANALYSIS 

In this section we provide an analysis of the X-ray, opti- 
cal, and radio light curves and spectra of the afterglow of 
GRB 050820A. We have divided the X-ray and optical light 
curves into segments (Phases 1^4) based on noticeable tempo- 
ral breaks. We then investigated each segment independently, 

20 The National Radio Astronomy Observatory is a facility of the Na- 
tional Science Foundation operated under cooperative agreement by Asso- 
ciated Universities, Inc. 
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FIG. 3.— fi-band imaging of the field of GRB 050820A. Left: P60 R c -band discovery image of the afterglow of GRB 050820A. The BAT (2'-radius white 
circle) and XRT (7"-radius black circle) error circles are labeled. The afterglow is identified with the two black tick marks. Center: The Second-Generation 
Digitized Sky Survey image of the identical field. The afterglow is not visible in this reference image. Right: Hubble Space Telescope F6 25W im age of the 
afterglow (indicated again with the two black tick marks). The two nearby objects complicating the ground-based photometry are visible (see 312.3. ll for details). 
All images are oriented with North up and East to the left. 



TABLE 4 

Radio Observations of GRB 050820A 



Observation Date 


fBAT 


Frequency 


Flux Density 3 


(2005 UT) 


(days) 


(GHz) 




Aug 20.39 


0.116 


4.86 


< 102 


Aug 20.39 


0.116 


8.46 


110±40 


Aug 20.39 


0.116 


22.5 


< 186 


Aug 21.20 


0.93 


8.46 


634 ±62 


Aug 22.42 


2.15 


4.86 


256 ±78 


Aug 22.42 


2.15 


8.46 


419 ±50 


Aug 22.42 


2.15 


22.5 


<216 


Aug 24.38 


4.11 


4.86 


171 ±47 


Aug 24.38 


4.11 


8.46 


74 ±36 


Aug 25.32 


5.05 


8.46 


< 114 


Aug 26.40 


6.13 


8.46 


< 120 


Aug 28.37 


8.10 


8.46 


166 ±45 


Sep 1.33 


12.06 


8.46 


89 ±39 


Sep 4.18 


14.91 


8.46 


106 ±33 


Sep 15.20 


25.93 


8.46 


76 ±30 


Oct 20.19 


60.92 


8.46 


< 70 



"Errors quoted for detections are at the lcr level. Upper 
limits are reported as 2<r RMS noise. 



fitting the light curve and spectra to power-law decay indices 
of the form: F v oc t~ a v~$ . The lack of a bright afterglow 
makes such analysis impossible in the radio. 

The X-ray and optical light curves, with temporal divisions 
marked as dashed vertical lines, are shown in Figure|4] Phase 
1 begins with the BAT trigger and ends with the resumption of 
X-ray observations at fb rea k,i = 4785 s. The X-ray and optical 
data behave differently in Phase 1, resulting in unique subdi- 
visions for the two band-passes. However, Phase 1 is the only 
epoch to show such divergent behavior. 

With the emergence of Swift from the SAA at /break, l, both 
the X-ray and optical light curves exhibit a relatively shallow 
decline. This characterizes Phase 2, which ends when the X- 
ray decline steepens at fbreak.2 = 8.7 x 10 4 s. The decay slope 
in Phase 3 is steeper than in Phase 2 in both band-passes. 

Phase 3 extends out to the last X-ray detection at fb re ak,3 = 
1.7 x 10 6 s. Between this time and the HST optical obser- 
vations, the optical decay must have significantly steepened. 
This last epoch, with only optical data, we define as Phase 4. 

The results of this power-law analysis are shown in Tables 



[3]and[5] Each bandpass is discussed in further detail below. 

3.1. X-ray Light Curve and Spectrum 

The X-ray light curve of GRB 050820A is shown in Fig- 
ure |4] (top panel). In Phase 1, we see two distinct behaviors: 
initially the X-ray light curve falls rapidly with a decay slope 
ctia.x = 2.2 ±0.3. This continues until tsAr = 217 s, and we de- 
fine this as Phase la x . The X-ray emission then rises rapidly 
(217 < ?bat < 257 s; Phase lb x ), after which Swift enters the 
SAA. It is clear from the correlation between the 7-ray and 
X-ray emission in this epoch that the two are related (see 
Fig . ^:) . A full study of the properties of Phase lb x is left 
to 

Following Swift's, emergence from the SAA, the light curve 
in Phase 2 shows evidence for a shallower epoch of decline. 
With the large gap in coverage, we cannot constrain when this 
transition occurs. We therefore define fb re ak,i to coincide with 
the resumption of XRT observations at ?bat = 4785 s. A sim- 
il ar br eak in the optical light curve is also seen near this time 

(€L3- 

The X-ray data after f break. 1 are not well fit by a single 
power-law decay (xj? = 3.7; 31 d.o.f.), due mostly to a steep- 
ening of the decay at ?bat ~ 10 5 s. Fitting a broken power- 
law model to this data, we find an acceptable fit with fbreak,2 = 
(8.7 ±2.4) x 10 4 s (x r 2 = 1.19; 25 d.o.f.). The resulting decay 
index before the break (Phase 2) is a.2. x = 0.93 ± 0.03. For 
Phase 3, we find a 3 , x = 1.25 ± 0.07. 

The X-ray spectral index in Phase la x is relatively steep: 
Aa,x = 0.90 ± 0.09. The spectrum hardens significantly in 
Phase lb x (/3ib, x = -0.10 ±0.03), further justifying our deci- 
sion to split Phase 1 into two separate X-ray segments. In 
Phase 2, the spectrum softens again, to (32.x = 1.20 ±0.04. 
There are too few X-ray counts in Phase 3 to meaningfully 
constrain the spectrum. 

The results of our analysis of the X-ray data set are summa- 
rized in Table 

3.2. Optical Light Curve and Spectrum 

The optical light curve from GRB 050820A is shown in Fig- 
ure|4](bottom panel). Phases 1 (0 < f BA T < 4785), 2 (4785 < 
f BA T < 8.7 x 10 4 ), and 3 (8.7 x 10 4 < t BAT < 1.7 x 10 6 ) of the 
optical light curve have already been defined in terms of the 
X-ray decay. However, unlike the X-ray, the earliest optical 
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TABLE 5 

Optical Afterglow Spectral and Temporal Fits 



Phase t«g t** a /3 X, 2 /d.o.f. 
(s) (s) 

la Tbat 626 -0.35 ±0.02 ... 3.00/ 13 

lb 626 4.8 xlO 3 0.97 ±0.01 0.57 ±0.06 1.53/ 14 

2 4.8xl0 3 8.7 X 10 4 0.78 ±0.01 0.77 ±0.08 5.7/31 

3 8.7 xlO 4 1.7xl0 6 0.99±0.06 ... 1.43/18 

4 1.7xl0 6 3.2 x 10 s >2.1 



NOTE. — We have fit the optical data to a power-law model of the form 
F v oc r a v~P where possible. In some Phases (la and 3), we have limited 
spectral coverage and could not meaningfully constrain (3. Thus we have only 
fit for the temporal decay index a. In Phase 4, we can only place an upper limit 
on the decay slope. All errors quoted are 90% confidence limits. 



observations indicate the afterglow was getting brighter with 
time (iCenko & Foxll2005l iWren et al] l2005). This rise con- 
tinues until a peak at ?bat ~ 600 s, marking the end of Phase 
la opt . After the peak, the optical light curve in all four P60 
filters decays steadily with aib, pt = 0.97 ±0.01 until fbreak.i- 
We note that for this and all subsequent phases, we have con- 
strained a to be identical in all opti cal fi lters. A more thor- 
ough discussion of Phase 1 is left to 34.21 

Much like the X-ray, the optical decay in Phases 2 and 
3 is poorly fit by a single power-law (Xr = 13-7; 49 d.o.f.). 
In Phase 2, the optical decay noticeably flattens (a2,opt = 
0.78 ±0.01). A much higher degree of variability is seen in 
the different filters, resulting in a poor fit statistic. After fbreak.2 
the decay in Phase 3 is again steeper and more uniform, with 
a 3 ,o P t = 0.99 ±0.06. 

It is clear that if we extrapolate the decay from Phase 3 out 
to the HST observations, the late-time flux is greatly overes- 
timated. We conclude therefore a break has occurred in the 
light curve sometime after fbreak,3> an d thus we define Phase 4 
to span 1.7 x 10 6 < f B AT < 3.2 x 10 6 . We estimate the tempo- 
ral decay in Phase 4 to be a4. opt > 2.1. 

Due to the limited spectral coverage of our observations, we 
are unable to provide meaningful constraints on the spectral 
index j3 in Phases la opt , 3, and 4. For the remaining epochs, 
we have excluded the U- and Z?-band data from our spectral 
fits, as these are expected to lie below the Ly-a absorption 
edge at this redshift. We attempted to solve for the host galaxy 
reddening (Ay(host)) using extinction laws fo r the Milky 
Way, Large and Small Magellanic Clouds from Pei ( 1992). 
In all cases, we find a host extinction consistent with zero. 
While this result is inconsistent with th e column density de - 
rived from high-resolution spectroscopy tLedoux et al.l2005h . 
it is in agreement with the low host extinction values seen 
in alm ost all well-sampled pie-Swift afterglows ( Ka nn et all 
2006). 

Ignoring host reddening, we find /3ib |0 pt = 0.57 ± 0.06 and 
Pi.opt = 0.77 ±0.08. While the optical spectrum appears to 
have steepened in Phase 3, the poor fit quality of this phase 
precludes any firm conclusions from being drawn. 

The results of our analysis of the optical data set are shown 
in TableU] 

3.3. Radio Light Curve 

The radio emission rises to a peak sometime around 1 d 
after the burst (see Figure [5] and Table |4}. The radio spec- 
trum at early times is quite chaotic, transitioning through a 



peak around 8.5 GHz at ?bat = 2.15 d to optically thin 4.11 
d after the burst. We note that some of the v ariation at 
early times may be due to inter-stellar scintillation (Goodman 
1997) as has been in seen in many other radio afterglows (see 
e.g. lFrail et al.ll997l) . 

The late-time (t > 7 d) radio data show no sign of any af- 
terglow brighter than 200 /iJy. This is in marked contrast to 
the bright optical and X-ray afterglows, and is on e of t he most 
difficult aspects of the afterglow to account for ( 34.4> . 

4. DISCUSSION 

In this section, we use the results from our previous 
analysis to try and explain the broadband emission from 
GRB 050820A in the context of the standard fireball model 
(see e.g. iPirad 120051 for a review). In this model, highly- 
relativistic ejecta are emitted by a central engine as shells with 
varying Lorentz factors. Shock fronts formed between these 
shells produce the high-energy prompt emission; these colli- 
sions are known as internal shocks. Emission from internal 
shocks is highly non-thermal and results in an (empirically 
determined) power-law spectrum with a cut-off exponential: 
dN/dE oc E a e- E l Eo (Ba nd et alJll993l) . As the relativistic 
shells expand and slow down, they eventually encounter the 
circum-burst medium. Again a collision-less shock front 
forms, accelerating electrons to a power-law distribution of 
Lorentz factors with exponent p and minimum Lorentz factor 
7 m . It is assumed that a constant fraction of the total energy 
density is partitioned to the electrons (e e ) and the magnetic 
field (cb)- These accelerated electrons then emit synchrotron 
radiation, powering the long-lived X-ray, optical, and radio 
afterglow (the forward, external shock). Additional emission 
can be generated from shock heating of the ejecta (the reverse 
shock), leading to a rapidly-decaying flare in the optical and 
radio light curve (Sari & Piran 1999b). 

The observed afterglow spectrum depends on the relative 
ordering of the three critical frequencies: v & , the frequency 
where self-absorption becomes important, v m , the character- 
istic frequency of the emission, and v c , the frequency above 
which electrons are able to cool efficiently through radia- 
tion. Typical afterglow observations occur when v a < v m < v c 
(i.e. slow cooling ), resulting in the following spectral indices 
ISariet all 19981) : 

F v oc v 1 ; v < v a (1) 

; v a < v < v m (2) 
V -{ P -\)I2 - jVm<v<Vc (3) 



8 



Cenko et al. 




Time Since Swift Trigger (s) Time Since Swift Tngger (a) 



FIG. 4. — X-ray (top) and optical (bottom) light curves of GRB 050820A. Both band-passes have been divided into four segments (Phases l^t), each shown 
with a vertical dashed line. The unique subdivision of Phase 1 is shown as a dotted line in both plots. Top: 2-10 keV X-ray fluxes were converted to flux densities 
at 5 keV using the average spectral slope for each Phase. It was assumed the spectrum remained constant from Phase 2 to Phase 3. The best-fit X-ray temporal 
decay is shown with a solid line. Bottom: Magnitudes were converted to flux densities using zero-points from Fukug ita et alj 119951) . The best-fit i?c" Dan d 
temporal decay is shown with a solid line. For full details of our analysis, see il3land Tablesl3landl5l 



v-pI 2 \v>v c (4) 

The light curve produced by such emission depends on 
the radial profile of the circum-burst medium into which the 
shock is expanding. The simplest circum-burst medium to 
consider is one in which the density is constant (p oc r°). This 
scenario is also referred to as an inter-ste llar medium (ISM). 
In this case, the flux density will scale as JSari et all!9 98): 

F v oc t '/ 2 ; v < v a (5) 
t '/ 2 ; v a < v < u m (6) 
1 3(1 " p)/4 ; v m < v < v e (7) 

f (2-3p)/4 . v > Vc (8) 

Alternatively, if we eliminate p from the above equations, we 
find a characteristic relation between a and in each sp ectral 
regime known as a "closure relation' ' dPrice et alJ2002k : 

a= f W<v A (9) 

f W,<v<v m (10) 

^ ; v m < v < u c (11) 

^\v>v, (12) 

The long-soft class of GRBs, however, is thought to arise 
from the deaths of massive stars, as they colla pse to form 
black-holes: the so-called "collapsar" model dWooslevll993l) . 
In the late stages of evolution, massive Wolf-Rayet stars are 
stripped of their outer envelopes in a wind, leaving behind a 
signature p oc r" 2 density profile that should be discernible in 
the afterglow light curve. The an alogous temporal decay in- 
dices for a wind-like medium are (Chevalier & Li 2000): 

F v cx t 1 ; v < v a (13) 

f° ; v a < v < v m (14) 

t (l -' ip) l A ;v m <v <u c (15) 

t (2-3p)/4 . v > Vc (16) 

The derived closure relations are: 

a= § ;v<v. (17) 

^L; V&<V <V m (18) 
^L; Vm<V <V, (19) 

; v > v, (20) 



The abo ve tem poral decay indices and closure relations 
(Equations I5I20> are only valid for spherically symmetric 
emission. GRBs, however, are thought to be beamed events. 
At early times, observers only notice emission from a nar- 
row cone (opening angle 9 ~ T" 1 , where T is the Lorentz fac- 
tor of the e xpanding shock) du e to relativistic beaming (see 
e.g. Rvbicki & Liahtman I1979I) . As the shock slows, how- 
ever, lateral spreading of the jet becomes important, and the 
observe r eventually notices "missing" emission from wider 
angles (lRhoadslll999t ISari et al1H99 9). This hydrodynamic 
transition manifests itself as an achromatic steepening in the 
afterglow light curve, with an expected post-jet-break decay 
proportional to t~ p . 

With the above formulation in hand, we now wish to under- 
stand the physical implications of our previous analysis. 

4. 1 . Early ~f-ray Emission 

The most striking feature of the 7-ray light curve of 
GRB 050820A is the large gap between the initial pulse that 
triggered the Swift BAT (Peak A in Fig. [TJ?) and the bulk 
of the high-energy emission (?bat > 200 s). The natural 
question arises as to whether this "precursor" 21 results from 
the same physical mechanism as the bulk of the high en- 
ergy emission. Many models predict a high-energy com- 
ponent distinct from the prompt GRB at early times. Pos- 
sible mechanisms include the transition from an optically 
thick to an optically thin environment in the fireball itself 
JPaczvnskilll98fl: iLvutikov & Usovll2000t iMeszaros & Reesl 
2000; Lvutikov & Blandford 2003), or the interaction of a jet 
with a progenitor, presumably a collapsing Wolf-Rayet star 
iRamirez-Ruiz et all2002t lWaxman & Meszaros 200j). 

We can securely rule out both of these models for the pre- 
cursor of GRB 050820A, as both predict a thermal spectrum. 
Fitting a thermal model to the BAT precursor spectrum results 
in a fit statistic of Xr =3.4 (75 d.o.f.), while a non-thermal 
power-law model provides an excellent fit (% 2 =1.07 for 75 
d.o.f.; T= 1.74 ±0.08). 

While the precursor may be non-thermal, it is noticeably 
softer than the majority of the remaining prompt emission 

21 Here we define a precursor as an event that is well separated from and 
contains only a small fraction of the total high-energy emission. Unlike some 
other authors, our definition is independent of the mechanism behind the 
emission. Peak A in GRB 050820A is then clearly a precursor. 
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FIG. 5.— Radio Afterglow of GRB 050820a. The early rise in the ra- 
dio light curve at (bat ~ 1 d is most easily understood as a reverse shock 
caused by late-time energy injection, as seen in both the optical and X-ray 
light curves. The most striking feature of the radio l ight curve, however, is 
the lack of a bright radio afterglow at late times (see j|4.41 . 
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FIG . 6 . — Joint BAT/XRT Spectrum of the Main Pulse of Prompt Emission 
(217 < (bat < 257 s). While there is no region of direct overlap, the XRT 
data is clearly the low-energy tail of the prompt BAT emission, forming one 
continuous spectrum. The best fit spectrum (T = 0.94) is shown as a solid 
line. Both the BAT and XRT data have been binned for plotting purposes. 



(see Figs. |2jl and|2^). A search for precursors in a sample of 
long, bright BATSE bursts revealed such a soft, non-thermal 
component in a sizable fraction (20-25%) of these events 
l lLazzatill2005t) . Furthermore, two of the lo ngest, brightest 
Swift bursts observed to date, GRB 041219 jVestrand et al l 
| 2005[lMcBreen et al.ll2006l) and GRB 060124 ( iRomano et aLl 
2006) show a faint, soft precursor followed by a large time lag 
(570 s in the case of GRB 060124). 

These soft precursors are inconsistent with the main prompt 
emission in most GRBs, which exhibit s a hard-to-soft evo- 
lution in the 7-ray spectrum (Ford et al. 1995; Frontera et al. 
2000). The 7-ray light curve of GRB 050820A conforms to 
this trend only if we ignore the precursor. Furthermore, it is 
difficult to conceive of a scenario by which internal shocks 
can generate such long periods of quiescence in a sustained 
outflow. The large time lag, soft nature, and repeated oc- 
currence of these precursors hint that they are in fact due to 
a different emission mechanism than the internal dissipation 
thought to power the bulk of the high-energy emission. How- 
ever, we cannot state this conclusively, as would be the case if 
the precursors were thermal. 

If we assume a different emission mechanism, the prompt 
emission did not begin until 222 s after the Swift trigger. This 
seemingly small discrepancy in defining 7b affects the cal- 
culated tem pora l decay indices, particularly during the early 
afterglow ( q4.4> . For all temporal decay indices calculated in 
this work, we consider 7b to coincide with the beginning of 
the bulk of the high-energy emission (i.e. 06:38:35 UT on 20 
August 2005). 

Finally, we consider the early X-ray emission. The tempo- 
ral decay slope at early times {a\ d .x = 2.2) is too steep to be at- 
tributed to a standard forward shock afterglow. The most pop- 
ular explanation for the rapid decline of early X-ray emission 
in Swift GRBs is "high-latitude emission": prompt emission 
from large angles (0 > T~ l ) that, due to relativistic beaming 
effects, reaches the o bserver at late times (Af ~ (1 +z)R0 2 /2c, 
iKumar & Panait escu 2000). However, this results in a well- 
defined relationship between the spectral and temporal indices 
(a = (3 + 2) which is inconsistent with the observed values for 
G RB 050820A 

Zhang et al. (2006) discuss possible mechanisms that could 
cause the early-time decay slope to deviate from this behavior. 
The most realistic possibility is if the X-rays were below the 
cooling frequency at this very early epoch. Then the closure 



relatio n would take the form a « 1 + 3/3/2 = 2.4 iSari & Piranl 
1999a), in good agreement with the observed value. 

4.2. Contemporaneous X-ray and Optical Emission 

Given the long duration and bright fluence, GRB 050820A 
provides a rare opportunity to study contemporaneous emis- 
sion in the optical, X-ray, and 7-ray band-passes. In Figure 
[fl we show the early time (?bat ^ 800 s) emission in X-rays 
(c) and optical (d) overlaid onto the high-energy light curve 
ofGRB050820A. 

A look at the X-ray data in Figure^ shows a strong correla- 
tion between the X-ray and 7-ray light curves. The X-ray light 
curve, previously in the midst of a decline, abruptly jumps in 
sync with the high energy emission at tsAT ~ 222 s (Phase 
lbx). In addition to temporal similarities, the X-ray photon 
index at this epoch (Txrt = 0.90 ± 0.03) is much harder than 
at any other epoch in the X-ray light curve, and similar to that 
derived from the BAT (Teat = 1 .07 ± 0.06). Thus motivated, 
we have performed a joint fit of the BAT and XRT spectra at 
this epoch. Unfortunately Konus-Wind had yet to trigger, and 
so no high-energy multi-channel spectra are available from 
that instrument. We find both band-passes well-fit by a single 
power-law with index T = 0.94 ±0.03 (x r 2 = 1-3; 391 d.o.f.). 
The resulting un-folded spectrum is shown in Figure [6] We 
conclude the X-rays in Phase lbx are generated by the same 
mechanism as the prompt emission. 

It is clear from Figure ^1 that, unlike in the X-ray band, 
there is no strong correlation between optical and 7-ray flux 
from GRB 050820A. Radical spectral evolution would be re- 
quired in the optical to explain both band-passes as arising 
from the same emission mechanism. We consider this sce- 
nario highly unlikely and conclude that, at the very least, the 
dominant contribution to the optical emission in Phase 1 has 
a different origin than the prompt emission. 

We next consider if our optical observations in Phase 1 can 
be explained solely in terms of the standard afterglow formu- 
latio n. We have attem pted to fit both a simple broken power- 
law ISariet alJl!998h as well as an analytic solution for the 
flux density near the optical peak ( Gran ot & Saril2002l) to our 
^c-, ^c~, and z'-band early-time data. The resulting fit quality 
is quite poor (\ 2 = 3.0; 14 d.o.f.) with the dominant contribu- 
tion coming from the data in Phase la opt (i.e. before the peak). 
In spite of the poor fit, we have included the results for Phase 
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TABLE 6 

Joint 7-ray / Optical Early-Time Data 



Interval ID tgjgi Duration 7-Ray Flux" a b F p b x? / d.o.f. Optical Filter Optical Flux Density 11 Co 7 c 

(s) (s) (1(T 8 ergcnr 2 s" 1 ) (keV) (mjy) 

1 80.949 85.376 < 1.3 ... ... ... V 0.20 ±0.03 < 10.2 

2 257.839 8.448 96.2+' 3 8 g 1.26±0.14 510+211 0.83/ 62 R c 2.17±0.08 12.5 

3 297.775 59.904 < 3.2 ... ..." ... / c 2.64 ±0.04 < 8.8 

4 389.167 57.344 113 ^om 269+^SJ 7 0.78 /62 z' 5.77±0.57 10.0 

5 487.471 57.344 9.35^ 1.96±0.18 ... 1.01/58 R c 4.28±0.08 9.2 

6 602.159 49.152 < 2.3 ... ... ... / c 4.14±0.05 < 7.9 

7 684.079 57.344 <4.1 ... ... ... z' 4.20 ±0.49 < 8.7 



a Errors quoted are at the 1 a level. 

b Spectral fits of the form dN/dE <x E~ a exp~' 2 ~ a ' £ / £ P were performed for the case of Intervals 2 and 4. For interval 5, the highest-energy data were 
not of sufficient quality to estimate F p . Instead a power-law fit (dN/dE oc E~ a ) was used. Errors quoted are 90% confidence limits. 
c 7-ray-to-optical color index: Co 7 = — 2.51og(F(opt)/F(7)). 



la opt in Table[5]for reference. 

This result is not unexpected, as IVestrand et alJ (2006) 
have shown that contemporaneous optical imaging of 
GRB 050820A with the RAPTOR telescope can be well-fit as 
the sum of two independent components: one representing 
the forward shock and another proportional to the high-energy 
prompt emission. We attempted an analogous fit with the P60 
and Konus-Wind data set. While a better fit statistic ensues, 
we still do not find an acceptable result = 2.2; 14 d.o.f.). 
We conclude that the relatively sparse time sampling of our 
observations, coupled with the frequent filter changes, make 
it impossible to verify this result. 

Independent of any correlation between the prompt optical 
and 7-ray emission, we note that the decay after ?bat = 600 s is 
dominated by the forward shock. Unlike the bright, early-time 
emission seen from GRB 990123 JAkerlof et all 119991) . we 
se e no evidence for rapidly-decaying (curs = (27p + 7)/35 « 
2; lKobavashii r2000) reverse-shock emission from an optical 
flare. 

Finally, using the combination of optical and 7-ray data, 
we consider the broadband spectral energy distribution (SED) 
of GRB 050820A at early times. For each of the contem- 
poraneous optical observations, we have extracted fluxes and 
(where possible) spectra from the corresponding Konus-Wind 
observations (F igure|7land Table|6j. Following the method of 
IVestrand etatl i2005). we have calculated the 7-ray-to-optical 
color index, Co 7 =-2.5 log(F(opt) /F(j)),oi lower limits, for 
each interval. The ratio varies significantly over the course of 
our observations. The value of Cq 1 = 12.5 in Interval 2 is 
consistent with that seen from GRB 041219a, while later in- 
tervals are even brighter in the optical. In fact, the optical-to- 
7-ray flux ratio in Interval 5 is over 240 t imes larger than that 
observed for GRB 050401 llRvkoffetal.1120051) . Evidently a 
large diversity exists in the broadband SEDs of GRBs at early 
times. 

4.3. Late-time Energy Injection 

The majority of XRT light curves observed to date have 
exhibited a period of shallow decline (0.2 < a < 0.8) 
that is inconsistent with the standard afterglow formula- 
tion (iNousek etalJ l2006). Two models have been invoked 
to explain this phase, both of which involve injecting en- 
ergy _jnto_the_forw<ffd shock at late times (f >> ?grb; see 
e.g. Zhang etaA 120061 for a review). In the first, the cen- 
tral engine is active for long time periods, t » ?grb- The 



late-time emission of highly relativistic material injects ad- 
dition al energy into the forward shock, flattening the decay 
slope dKatz & Piranlll997llRees & Meszarosll2000l) . Alterna- 
tively, towards the end of the 7-ray emission, the central en- 
gine may inject material with a smooth distribution of (de- 
creasing) Lorentz factors. Slower moving material will catch 
up with the forward shock when it has swept up enough 
circum-burst material, resultin g in a smooth injection of en- 
ergy at late times (iRees & Meszarosll 19981 ISari & Meszaros 
2000). While both models explain the flattening of the XRT 
light curves, they provide different constraints on the nature 
of the central engine. 

In the first (long-lived central engine) model, the central 
engine's luminosity, L(t), is characterized as 



L(t): 



(21) 



This results in the following spectral and temporal power-law 
indices for a constant-density medium: 

5^-8 



V < V xx 
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For a wind-like medium, the analogous results are 
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FIG. 7. — Early-time Broadband Spectral Energy Distribution of GRB 050820a. Konus-Wind spectral data for the given intervals (see Tablelfjlfor definitions 
of the intervals) are shown (crosses) alongside the corresponding optical observations (circles). 2-sigma upper-limits in the high-energy spectra are plotted as 
triangles. The best-fit model to the Konus-Wind spectrum is shown as a solid line, and the dashed lines show the 90% confidence intervals for the spectral 
fits. The ratio of optical to 7-ray flux varies significantly between the three intervals. Left: Interval 2. This interval covers only a small fraction of the time 
of the corresponding P60 image because Konus-Wind only triggered on GRB 050820a (and hence began collecting multi-channel spectra) in the middle of this 
image. The optical Rf-band data point lies below the extrapolation of the 7-ray spectrum. Center: Interval 4. Here the Konus-Wind and P60 intervals are nearly 
simultaneous. Unlike the other intervals, the P60 z'-band point lies above the predicted value and within the 90% confidence interval of the extrapolation of the 
high-energy spectrum. Right: Interval 5. Here the high-energy extrapolation greatly overestimates the optical flux. However, the 7-ray flux in this interval is 
quite low, and in particular the low number of high-energy photons makes it difficult to constrain a cut-off power-law spectrum. In fact this interval was best fit 
with a simple power-law spectrum. 



q-2 | (2 + q)(3 



(27) 



The refreshed shock scenario is parameterized in terms of 
the amount of mass ejected with Lorentz factor greater than 



7: 



M( > 7) oc 7 



(28) 



For the circum-burst profiles considered here, we can define 
a new variable, q, such that w e reprod uce identical aftergl ow 
behavior to that of Equations 1221241 or Equations 1251271 by 
simply substituting q for q. q is relate d to the mass eject ion 
parameter s by the following equations (Zhang et al. 2006): 



(ISM) 



3+5 



(Wind) 



(29) 
(30) 



While the X-ray decay in Phase 2 is not as fiat as that seen in 
other Swift bursts, the temporal and spectral decay indices are 
nonetheless inconsistent with the standard afterglow model 
for z^(X) > v c (Equations 1121 and I20>. Furthermore, the op- 
tical light curve shows a flattening during Phase 2, and is in- 
consistent with both the closure relations in either medium for 
y m < ^(opt) < v c (Eauationsll lland!19ll. We conclude we are 
therefore seeing a milder version of the energy injection phase 
present in many Swift XRT afterglows. 

For the X-ray data in Phase 2, we find an acceptable fit 
for the energy injection models only if ^(X) > v c . This cor- 
responds to values of qx = 0.66 ± 0.08 (sx,ism = 2.0 ± 0.3, 
sx.wind = 3.1 ±0.7) and p x = 2.4 ±0.2. The optical data 
in Phase 2 are best fit with a constant-density medium and 
v m < v(opt) < u c : q ovt = 0.73 ± 0.09 (si S M, op t = 1.8 ±0.3) and 
p pt = 2.5 ±0.2. Both the X-ray and the optical fall in the spec- 
tral regime we would expect, providing further confidence in 
this interpretation. 

A prediction of the energy injection hypothesis is a bright 
reverse shock at early t imes most easily visible in the radio 
ISari & Meszarosl2000l) . A reverse shock nicely explains the 
rapid decline in flux at 8.5 GHz from 1-4 d after the burst. 
Furthermore, the transition from a spectrum peaked around 
8.5 GHz at f BAT = 2.15 d to an optically thin radio spectrum at 
?bat = 4.11 d can be understood as the reverse shock peak fre- 
quency, i^ s , passing through the radio. Since » f^/7 2 , 



this should occur well before the forward shock peak fre- 
quency reaches the radio bands. 

Distinguishing between the two theories to explain the en- 
ergy injection is quite difficult, as both models can be iden- 
tically parameterized. Progress in this area would require a 
large sample of bursts with detailed contemporaneous X-ray 
and optical light curves. If the refreshed shocks are due to 
continued engine activity, they should be correlated with the 
bright X-ray flares seen in some XRT afterglows. On the other 
hand, if the flat decay is caused by slow-moving ejecta, this 
behavior should be more uniform from burst to burst. Such an 
analysis is beyond the scope of this work. 

4.4. Burst Environment and Progenitor Models 

We now turn our attention to the issue of the circum-burst 
medium. As discussed earlier, the radial profile of the burst 
environment affect s the tempor al decay below the cooling fre- 
quency (EQuations l5l7irT3ll5>. In particular, the closure re- 
lationships (Equations 1911 II 1171191 are sufficiently different 
that we should be able to distinguish between the competing 
models for a well-sampled event like GRB 050820A. 

Firs t we examine the X-ray data. As discussed previously 
( 34.31 1. the X-ray observations in Phase 2 require invoking 
mild energy injection to explain the shallower-than-expected 
decay for ^(X) > v c . If we assume the X-ray spectral index 
does not change from Phase 2 to Phase 3, then we find «3,x 
and (3i.x satisfy the standard afterglow closure relation for 
vx > v z (Eauationsll2landl20>. The corresponding values for 
the electron index are: p a , x = 2.3 ± 0.1, pp ix = 2.40 ±0.08. 

Unlike the X-ray observations, the optical bands typically 
probe frequencies below the cooling frequency, where the clo- 
sure relations are different for different circum-bur st m edia 
(Equations 1111 and 1 1 91. We have shown already in 34.31 that 
the optical data in Phase 2 are better fit by a constant-density 
medium. We find that a constant-density medium is favored 
in the optical in Phase lb op , and Phase 3 as well. The only 
closure relation satisfied in Phase lb opt is for an ISM with 
v m < Kopt) < v z (Eauationll 1>. The resulting p-values are 
p Qlbopl = 2.29 ±0.02, Pf3lb v = 2.1 ±0.1. We note that had we 
equated the BAT trigger time, That, with the onset of the burst 
(To), the temporal slope in Phase lb opt would not have been 
consistent with any closure relation. 
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In Phase 3, we cannot meaningfully constrain the optical 
spectral slope. However, using the X-ray-to-optical spectral 
slope in this phase, /3ox,3 ~ 0.8, we conclude the optical data 
in this segment still fall below the cooling frequency. Based 
solely on the temporal decline then, we can rule out a wind- 
like medium in this phase, The corresponding electron index 
(p w 1 .7 ± 0. 1) would result in a divergent total energy. While 
this possibility has been addressed with more c omplicated 
electron energy distributions (see e.g. Dai & Cheng 2001), we 
consider this possibility unlikely. 

Taken together, the X-ray and optical data provide a con- 
sistent picture of the forward shock expanding into a constant 
density medium. The late-time (/bat > 1 week) radio observa- 
tions, however, are inconsistent with this interpretation. For a 
constant-density medium, the peak flux density, F v , max , should 
remain constant in time. This would predict, if we have cor- 
rectly interpreted the optical peak as the forward shock ( §I4.2> . 
a similar peak (F vrmx w 5 mJy) in the radio at £bat ~ 7 days 
(v m oc f" 3 / 2 ). This is well above the VLA detection limit at 
this epoch, yet we only measure F v ~ 100 /iJy. While the 
energy injection phase will delay the arrival of v m in the ra- 
dio (fm Jec oc f" 3 / 2 f 3 ( s-I )/ 2 ( 7+J )), our radio limits extend out to 
two months after the burst. It would be very difficult, if not 
impossible, to delay the peak this long. Furthermore, dur- 
ing the energy injection phase, the peak flux increases with 
time (/v"max oc f 3 ( I_I )/( 7+I )). Thus we would expect to see rising 
emission earlier relative to the peak, counteracting the delay 
of the peak radio flux. 

One explanation for the lack of a bright, late-time radio af- 
terglow is an early jet break (t < 1 day), as was invoked for 
GRB 990123 (Kulk arni etalJ 1999). However, we find no ev- 
idence for a jet break in the optical or X-r ay light curves out 
to at least 17 days after the burst (see j!4.5> , 

Another possibility, invoked to explain the relatively low 
late-ti me radio flux fr om GRB 050904, is a high ambient den- 
sity JFrail et al J 120061) . In the case of GRB 050904, it was 
argued that the large density raised the self-absorption fre- 
quency, v a , above the radio observing bands. This greatly 
suppresses the radio flux, for the spectrum in this regime is 
proportional to v 2 (Eqn.Q. There is no evidence in the radio 
data for an optically thick spectrum, although spectral data is 
sparse at lat e times. Furthermore, broadband modeling of this 
event ( j34.5> rules out a high ambient density for typical val- 
ues of the micro-physical parameters e e and e^- We therefore 
consider this explanation unlikely. 

Alternatively, a natural explanation for the low radio flux at 
late times is a wind-like medium. In a wind-like medium, the 
forward shock peak flux density declines in time as F„ )max oc 
f -3/2 rp ne decreasing peak flux counteracts the rising syn- 
chrotron emission, suppressing any late-time radio data. This 
is of course inconsistent with our X-ray and optical data, 
which strongly favor a constant-density medium. One can 
imagine a scenario in which the environment near the burst 
(the regime sampled predominantly by the X-ray and optical 
data) is approximately constant in density, while the outer re- 
gions (sampled by the radio at later times) have a wind-like 
profile. However, without a physical justification for such a 
density profile, this remains little more than speculation. The 
lack of a bright radio afterglow remains a puzzling aspect of 
GRB 050820A. 

4.5. Geometry and Energetics 



Using t he hig h-energy fluence derived from Konus-Wind 
satellite ( 32.1.21 . we calculate the total isotropic energy re- 
lease in the prompt emission was £V i sn = 8.3* 2, i x 10 53 ergs 
(assuming a redsh ift of z = 2.615; IProchaska et al.l 120051 
iLedoux et aLll2005l) . This makes GRB 050820A one of the 
most energetic events (in terms of E^^) for which a redshift 
has been measured ( Amati 2006). 

However, only a fraction of the explosion energy is con- 
verted into prompt emission via internal dissipation. The rest 
remains in the kinetic energy of the outflow, powering the for- 
ward shock and hence the afterglow. We can estimate the ki- 
netic energy of the afterglow ( £kE.iso) by exam ining the the X- 
ray emission at £bat > 10 hr dFreedman & Wa xman 20011) . At 
this point, the X-rays should be above the cooling frequency. 
The flux density is then independent of ambient density and 
only weakly dependent on e#. A joint fit of the Phase 3 opti- 
cal and X-ray data (after the energy injection has stopped and 
the system has returned to adiabatic expansion) constrains the 
electron energy index: p = 2.34 ± 0.06. If we take typical val- 
ues for t e (0.1 -0.3) and e B (0.01-0.1) (lYost et all2003l) . we 
find that 15 < EicE,i S o,52 £ 100. 

For an accurate accounting of the total energy emitted by 
this event, however, we must determine the degree of colli- 
mation of the emission. We therefore examine all the tempo- 
ral breaks in the optical and X-ray light curves to determine 
which one (if any) shows an achromatic ste epening to the t~ p 
decay expected from a jet (Sa ri et aO 19991) . The only plausi- 
ble candidate is the transition from Phase 3 to 4 in the optical 
light curve. The steepening here is achromatic (i.e. it is seen in 
all three Hubble Space Telescope filters) and much too large to 
be explained solely by the cooling frequency passing through 
the optical bands (although this may have occured as well). 
Any contribution from an underlying host galaxy would only 
further steepen the decay in Phase 4. 

With only one observation, it is impossible to constrain the 
post-break decay index. Instead, we assume the post-break 
decay has a power-law index a = p w 2.34 (see above). We 
then find *j et< j = 18 ±2 d. This result is consistent with our 
X-ray observations, which put a lower limit on the jet break 
time of fj et .d > 17 d. 

We note that the jet break time we have inferred for 
GRB 050820A is extremely large. In the host galaxy refer- 
ence frame, the break occurs at ffe^ w 5, a factor of three 
larger than any jet break seen in the pre-Swift era (IZeh et alJ 
2006). In this respect, too, GRB 050820A is a strong outlier. 

To convert the jet break times to a range of opening angles, 
we use the relation (lSarietalJl999h : 

= 0.161 (-U (J?bu) (31) 

\l+Zj \ £ 7,iso,52/ 

Here ry 7 is the fraction of the total energy converted to prompt 
7-ray emission. The only remaining unknown in EquationOll 
is the ambient density, n. Using the ratio of the X-ray and 
optical data, as well as the canonical values of e e and €b, we 
find the density of the burst must be low: n < 1 cm" 3 . After- 
glow modeling of the late-time optical and X-ray data (Phase 
3, after any continued energy injection has ceas ed and the 
shock expands adiabatically) using the technique o flYost et alJ 
(2003) confirms this result: n < 0.1. 

Combining the above results, we find the opening angle is 
constrained to fall between between 6.8° < 9 < 9.3°, corre- 
sponding to a beaming factor /j, = 1 - cos 8 w 10~ 2 . While 
the opening angle is large for a long-soft burst, there are 
several comparable events in the pre-Swift sample (IZeh et alJ 
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FIG. 8. — Total energy release of GRB050820A. Left: Collimation-corrected energy release in the prompt emission (E~,) of a sample of cosmological 
GRBs, including GRB 050820 A. Right: Collimation-corrected blas t -wave energ y (£ke) for the same sa mple. (References: Bereer et al. 2004; Yost et all2003t 
Panaitescu & Kumar 2002; Berger et al. 2001; Chevalier et al. 2004; Soderbere et al. 2004; Bereer et al. 2003b). 
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2006). The total collimation-corrected energy emitted in 7- 
rays (£ 7 ) from GRB 050820A is therefore 7.5^ x 10 51 ergs. 
The corresponding value for the blast-wave energy (Eke) is 
5 2 + I ? x 10 S1 ergs . 

iGhirlanda et all (|2004a ) have demonstrated an empirical re- 
lation between E 1 and the peak energy of the prompt emis- 
sion spectrum in the GRB rest frame (E r p est ). GRB 050820A 
is more energetic than any of the 37 bursts considered in their 
sample (in terms of £ 7 ), and therefore proves an interesting 
test case for the so-called Ghirlanda relation. Using our calcu- 
lated value of E y , Ghirlanda predicts E" 

is marginally consistent with the actual value of E™ 1 
MeV derived from the Konus-Wind dataset. 

In Figure [8] we plot a histogram of £ 7 and Eke for the 
~ 15 long-soft cosmological bursts for which both quantities 
have been derived. We have not included the most nearby 
events (GRB 980425, GRB 031203, and GRB 060218) in our 
analysis, as these events released si g 
than the typical cosmological GRB (Soderberg et alJl20~O 
Soderberg et al. (2004) have shown that, with the exception 
of the most nearby events, the sum of E 7 and Eke is clus- 
tered around 2 x 10 51 ergs. GRB 050820A is clearly an over- 
energetic exception, an order of magnitude more energetic 
than this sample. In fact, it would require the direct con- 
version of « 10~ 2 M Q (with 100% efficiency) to release this 
much energy. 

Finally, it is important to consider how robust our estimates 
of £KE,isoi 0, and n are given that the standard afterglow model 
fails to explain the broadband behavior of GRB 050820A. We 
note that the opening angle is relatively insensitive to both 
variables (Ean.l31>: factors of order unity will be greatly re- 
duced by the ^ exponent. It is difficult to conceive of a long- 
soft GRB environment where the ambient density is less than 
10~ 2 , and high densities would only increase the opening an- 
gle and thus the energy release. 

5. CONCLUSION 

GRB 050820A joins a select sample of events with simul- 
taneous observations in the 7-ray and optical bands, and an 
even smaller group with contemporaneous X-ray observa- 
tions as well. Such events have led to fundamental advances 
in our understanding of GRBs, including the discovery of a 
revers e-shock optical flash from GRB 990123 ( Akerl oTet al.l 
1999) and possibly from GRB 050904 jBoer et all2006l) . The 



early-time optical emission from GRB 0412 19A also showed 
a bright optical flash, but the rise was correlated with an ac- 
companying peak of 7-ray e mission, su ggesting a common 
origin for the two components ( Vestrand et a lJ2005l) . 

The early behavior of GRB 050 820A is unlike either of 
these events. Vestrand et al. (2006) have shown the contem- 
poraneous optical emission is well described as the sum of 
two components: one proportional to the prompt 7-ray emis- 
sion and one smoothly varying forward shock term. While 
the 7-ray component is important for Tbat < 300 s, the op- 
tical peak at ?bat ~ 600 s is dominated by emission from 
the external shock region. Furthermore, the post-peak de- 
cay rate is inconsistent with reverse-shock emission. Instead 
we interpret this as the forward shock peak frequency pass- 
ing through the optical bands. Thi s is not unlike what was 
seen in the optical for GRB 060124 JRomano et all2006l) . al- 
though the time resolution in the prompt phase was much 
poorer than for this burst. The contemporaneo us optical light 
curve s of GRB 050319 (Oui mbv et afll2006l IWofniak et al] 
120051 and GRB 050401 ( Rvkoff et al 120051) did not show this 
peak phase, but extrapolations to late-times were consistent 
with the adiabatic expansion of a forward shock. Another dif- 
ferent behavior was seen in the early optical light curve of 
GRB 050801, which showed an extended plateau phase cor- 
related with the X-ray emission, hinting at continued energy 
injection from a cen tral engine refreshing the external shocks 
( Rvko ff etalJ2003l . 

The contemporaneous X-ray emission, on the other hand, is 
the low-energy tail of the p rompt emission. This behavior was 
also seen for GRB 060124 ( Roman o et alJ200 6) and has been 
hinted at in the rapid decline in e arly X-ray light c urves at- 
tributed to high-latitude emission (Liang et al. 2006), as well 
as the bright X-ray fl ares seen in many XRT light curves 
(Burr ows et al]l2005bl) . It is clear then, that, in marked con- 
trast to the X-ray emission, contemporaneous optical emis- 
sion exhibits a large diversity in behavior. Unfortunately the 
physical mechanism behind this diversity remains to be ex- 
plained. 

The issue of burst geometry is a particularly interesting one 
in the Swift era. The steep post-break decay slope, seen simul- 
taneously in multiple filters, makes GRB 050820A one of the 
most convincing examples of a beamed event in the Swift sam- 
ple. The X-ray afterglow, however, is too faint at late times 
to provide broadband confirmation. In fact, very few Swift 
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bursts, including those, like GRB 050408 JFolev et alJl2006l) . 
that have been followe d for months, show s igns of a jet break 
in the XRT light curve jNousek etall20 06). 

Typical jet breaks i n pie-Swift burs ts occurred on time- 
scales of several days iZeh et al] 120061) . iPerna et all fc003l) 
predicted that Swift would detect bursts with wider opening 
angles than previous missions due to the increased sensitiv- 
ity of the BAT. However, not a single Swift afterglow has 
shown a convincing jet-break-like transition in multiple band- 
pas ses (candidates includ e GRB 050525 A and GRB 050801; 
see Panaitescu et al. 2006 and references therein for a more 
thorough discussion). While it may be that most jet breaks, 
like GRB050820A, occur at late times, beyond the sensi- 
tivity of the XRT and most ground-based facilities, this is 
nonetheless troubling. On the one hand, few if any X-ray 
jet breaks were seen in pre-Swift bursts; all collimation an- 
gles were determined from optical (and sometimes radio) light 
curves. Conversely, given the large number of well-sampled 
XRT light curves, and the fact that such fundamental results 
for GRB cosmology as the Ghirlanda relation rest on our pic- 
ture of GRBs as a-spherical events, this is clearly a matter that 
merits further investigation. 

One consequence of the large opening angle associated 
with GRB 050820A is a correspondingly large burst and 
aft erglow energy. In fact , of all the bursts compiled in 
the [Ghirlanda etjjl] J2004bl) sample, GRB 050820A has the 
larges t prompt energy rele ase. And unlike GRB 990123 
( Panai tescu & K umar 2001), this large 7-ray energy was ac- 
companied by a correspondingly large kinetic energy im- 
parted to the afterglow. The only comparable event for which 
such energies could be determined was the high-redshift 
GRB 050904, which re leased a total en ergy of ~ 10 52 ergs 
(Tagliaferri et al. 20051 [Frail et al.l l2006). Given the large 7- 
ray fiuence, similar events should have been easily detected 
by both Swift and previous GRB missions. And given the 
bright optical afterglow and the late jet-break, such events are 
strongly favored for ground-based follow-up (i.e. redshift de- 
termination). The lack of a large sample of such events means 
they must be relatively rare in the universe. 

Like many other Swift GRBs, the X-ray light curve of 
GRB 050820A exhibits a phase of shallow decay incompat- 
ible w ith the standard forward shock model (Nous ek et alJ 
2006). GRB 050820A is relatively unique, however, in that 
this epoch is also well-sampled in the optical. The seem- 
ingly simultaneous breaks in the optical light curve bolster 
the commonly-held belief that this phase is caused by some 
form of refreshed shocks (Zhang et al. 2006). Coupled with 
the large gap between the precursor and the bulk of the prompt 



emission, the late-time energy injection poses fundamental 
challenges to any central engine model. 

Finally, we return to the question of the radio afterglow. 
Radio observations typically probe low-Lorentz factor ejecta 
(r ~ 2-3) at large distances from the central engine (r ~ 10 17 
cm). The forward shock peak frequency reaches the radio 
much later than the optical or X-ray bands. Thus, radio emis- 
sion is usually visible at later times than optical or X-ray emis- 
sion, and is well-suited to study afterglows when the emis- 
sion is isotropic (i.e. after the jet break) or even in some 
cases when the eje cta has slowed to Newtonian expansion 
(Ber reretalJ l2004). For this reason, late-time radio obser- 
vations are considered the most accurate method for model- 
independent calorimetry. For GRB 050820A, this paradigm 
has broken down. The burst had a bright optical and X-ray 
afterglow, but weak emission in the radio. It is hoped that fur- 
ther studies of such energetic GRBs in the Swift era will help 
to elucidate some of these puzzles. 
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Cenko et al. 



TABLE 7 

XRT Observations of GRB 050820A 



— : 

Mean Observation Date 


f BAT 


Duration 


Spectral Index (p) 


2-10 kev Mux 


(2005 UT) 


(s) 


(s) 




(10" ergs cm -2 s~') 


Antr 90 0fv3fv?S 

rtUg Z.U UU.JU.ilJ 


92.0 


10 


90 


14.9 ±1.3 


Aug 20 06:36:35 


102.0 


10.0 




12.2± 1.2 


Anp 20 06-36 45 


1 12.0 


10.0 




10.8 ±1.1 


AncT 90 06-36-55 


122.0 


10.0 




9.8 ± 1.1 


An? 20 06-37 05 


132.0 


10.0 




7.3 ± 0.9 


Aug 20 06:37:15 


142.0 


10.0 




8.1 ± 1.0 


Aug 20 06:37:25 


152.0 


10.0 




4.8 ±0.8 


Aug 20 06:37:35 


162.0 


10.0 




4.5 ±0.8 


Aug 20 06:37:45 


172.0 


10.0 




4.9 ±0.8 


Aug 20 06:37:55 


182.0 


10.0 




2.5 ±0.6 


Aug 20 06:38:05 


192.0 


10.0 




4.0 ±0.7 


Aug 20 06:38:15 


202.0 


10.0 




2.3 ±0.8 


Aug 20 06:38:25 


212.0 


10.0 




2.4 ± 0.7 


Aug 2U Uo.jo.jj 


222.0 


10.0 


-0. 10 


i /j.i ± /.y 


Alio ?n 


232 


10 






Aug 20 06:38:55 


242.0 


10.0 




631.3 ± 15.0 


Aug 20 06:39:05 


252.0 


10.0 




629.7 ± 26.5 


Aug 20 07:56:43 


4.910 X 10 J 


250.0 


1.20 


3.0 ±0.2 


Aug 20 08:00:53 


5.160 x 10 3 


250.0 




3.1 ±0.2 


Aug 20 08:05:03 


5.410 x 10 3 


250.0 




3.2±0.2 


Aug 20 08:09:13 


5.660 x 10 3 


250.0 




2.8 ±0.2 


Aug 20 08:13:23 


5.901 x 10 3 


250.0 




2.6 ±0.2 


Aug 20 11:41:43 


1.841 x 10 4 


2.5 X 10 3 




1.02 ±0.05 


Aug 20 12:23:23 


2.091 x 10 4 


2.5 X 10 3 




0.9 ±0.1 


Aug 20 13:05:03 


2.341 x 10 4 


2.5 X 10 3 




0.71 ±0.03 


Aug 20 14:28:23 


2.841 x 10 4 


2.5 X 10 3 




0.58 ±0.03 


Aug 20 15:10:03 


3.091 x 10 4 


2.5 X 10 3 




0.51 ±0.06 


Aug 20 15:51:43 


3.341 x 10 4 


2.5 X 10 3 




0.52 ±0.04 


Aug 20 16:33:23 


3.591 x 10 4 


2.5 X 10 3 




0.56 ±0.06 


Aug 20 17:15:03 


3.841 x 10 4 


2.5 X 10 3 




0.47 ± 0.07 


Aug 20 17:56:43 


4.091 x 10 4 


2.5 X 10 3 




0.49 ± 0.05 


Aug 20 19:20:03 


4.951 x 10 4 


2.5 x 10 3 




0.39 ±0.08 


Aug 21 15:48:01 


1.196 x 10 3 


1.0 x 10 4 




0.18 ±0.03 


Aug 21 18:34:41 


1.296 x 10 s 


1.0 x 10 4 




0.14±0.01 


Aug 21 21:21:21 


1.396 x 10 5 


1.0 x 10 4 




0.13 ±0.02 


Aug 22 00:08:01 


1.496 x 10 5 


1.0 x 10 4 




0.10 ±0.02 


Aug 22 02:54:41 


1.596 x 10 5 


1.0 x 10 4 




0.07 ±0.01 


Aug 23 14:30:22 


2.877 X 10 5 


1.0 x 10 4 




0.046 ±0.006 


Aug 23 17:17:02 


2.977 x 10 5 


1.0 x 10 4 




0.050 ±0.008 


Aug 24 18:02:47 


3.869 x 10 5 


3.5 x 10 4 




0.034 ±0.004 


Aug 25 03:46:07 


4.219 x 10 5 


3.5 x 10 4 




0.05 ±0.01 


Aug 27 19:38:46 


6.518 x 10 5 


1.25 x 10 5 




0.022 ±0.003 


Aug 29 06:22:06 


7.768 x 10 5 


3.5 x 10 4 




0.014 ±0.002 


Sep 04 15:45:23 


1.329 x 10 6 


1.0 x 10 5 




(6.8 ±2.5) x Kr 3 


Sep 05 19:32:03 


1.429 x 10 6 


1.0 x 10 5 




(4.2±1.5)x Kr 3 


Sep 06 23:18:43 


1.529 x 10 6 


1.0 x 10 5 




(8.8 ±3.3) x Kr 3 



NOTE. — The four Phases of the X-ray light curve are delineated by horizontal lines (see i!3.1l for 
further details). We assumed the spectral index was constant in each Phase to convert count rates to the 
flux values shown here. We also assumed the spectral index remained constant from Phase 2 to Phase 
3. All errors quoted are at the lcr level. 
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TABLE 8 

Optical Observations of GRB 050820A 



A/Tpnn Oh*;prvfitinn Oatp 

IVIGCIU L/USCl VaLLUIl l—'Cllt 


'BAT 




Filter 


pYnn<;nrp Timp 


A/T fi o n i ti i rl p a 


I? pfpi-pnpp 


(2005 UT) 


(s) 






(si 


CVeaa) 




Aug 20 08:14:47 


5.994 x 10 3 


UVOT 


U 


693.3 


18.12 ±0.13 




Aug 20 12:49:02 


2.245 x 10 4 


UVOT 


u 


899.8 


19.11 ±0.14 




Aug 20 16:01:58 


3.404 x 10 4 


UVOT 


u 


899.8 


19.40 ±0.14 


* 


Aug 20 19:14:54 


4.560 x 10 4 


UVOT 


u 


899.8 


19.66 ±0.15 


* 


Aug 20 22:59:51 


5.910 x 10 4 


UVOT 


u 


392.3 


19.70 ±0.18 


* 


Aug 20 09:36:06 


1.087 x 10 4 


UVOT 


B 


899.8 


18.79 ±0.16 


* 


Aug 20 13:04:10 


2.336 x 10 4 


UVOT 


B 


899.8 


19.24 ±0.15 


* 


Aug 20 16:17:05 


3.493 x 10 4 


UVOT 


B 


897.3 


19.74 ±0.16 


=1= 


Aug 20 19:30:02 


4.651 x 10 4 


UVOT 


B 


899.8 


19.87 ±0.16 


# 


Aug 21 23:30:41 


1.473 x 10 5 


RTT150 


B 


1800.0 


21.28 ±0.06 


1 


Aug 22 22:58:17 


2.318 x 10 5 


RTT150 


B 


7860.0 


22.05 ± 0.06 


2 


Aug 23 22:09:05 


3.153 x 10 5 


RTT150 


B 


5400.0 


22.38 ± 0.08 


2 


Aug 20 07:01:53 


1.620 x 10 3 


P60 


g 


120.0 


16.27 ± 0.04 


* 


Aug 20 07:12:35 


2.262 x 10 3 


P60 


g 


120.0 


16.65 ± 0.04 


=i= 


Aug 20 07:23:24 


2.911 x 10 3 


P60 


g 


120.0 


16.92 ± 0.05 


* 


Aug 20 08:04:47 


5.394 x 10 3 


P60 


g 


360.0 


17.59 ± 0.05 


* 


Aug 20 08:29:42 


6.889 x 10 3 


P60 


g 


360.0 


17.87 ± 0.06 


* 


Aug 20 08:54:50 


8.397 x 10 3 


P60 


g 


360.0 


18.11 ±0.06 


* 


Aug 20 09:33:18 


1.071 x 10 4 


P60 


g 


720.0 


18.29 ±0.06 




Aug 20 10:43:25 


1.491 x 10 4 


P60 


g 


720.0 


18.36 ±0.06 


* 


Aug 20 11:39:32 


1.829 x 10 4 


P60 


g 


720.0 


18.59 ± 0.09 


* 


Aug 20 06:36:58 


125.0 


UVOT 


V 


89.0 


18.18 ±0.20 


=1= 


Aug 20 07:52:58 


4.685 x 10 3 


UVOT 


V 


99.8 


17.20 ±0.13 


* 


Aug 20 10:04:28 


1.258 x 10 4 


UVOT 


V 


337.0 


18.16 ± 0.12 


* 


Aug 20 11:12:34 


1.666 x 10 4 


UVOT 


V 


899.8 


18.39 ± 0.09 


* 


Aug 20 14:25:30 


2.824 x 10 4 


UVOT 


V 


899.8 


18.86 ± 0.11 


* 


Aug 20 17:38:26 


3.981 x 10 4 


UVOT 


V 


899.8 


19.35 ± 0.12 


* 


Aug 20 21:03:50 


5.214 x 10 4 


UVOT 


V 


899.8 


19.14 ±0.11 


* 


Aug 21 16:09:34 


1.209 x 10 5 


UVOT 


V 


899.8 


20.24 ±0.19 




Aug 21 19:22:30 


1.325 x 10 5 


UVOT 


V 


899.8 


20.40 ± 0.22 




Aug 20 06:38:49 


236.0 


P60 


Rc 


60.0 


15.39 ± 0.04 


* 


Aug 20 06:43:31 


517.0 


P60 


Rc 


60.0 


14.65 ± 0.02 


* 


Aug 20 06:48:25 


812.0 


P60 


Rc 


60.0 


15.06 ± 0.03 


* 


Aug 20 06:53:59 


1.146 x 10 3 


P60 


Rc 


120.0 


15.42 ± 0.02 


* 


Aug 20 07:04:34 


1.781 x 10 3 


P60 


Rc 


120.0 


16.05 ± 0.03 


* 


Aug 20 07:15:19 


2.426 x 10 3 


P60 


Rc 


120.0 


16.39 ± 0.02 


* 


Aug 20 07:32:13 


3.440 x 10 3 


P60 


Rc 


360.0 


16.77 ± 0.02 


* 


Aug 20 08:12:58 


5.885 x 10 3 


P60 


Rc 


360.0 


17.32 ± 0.04 


* 


Aug 20 08:38:05 


7.392 x 10 3 


P60 


Rc 


360.0 


17.48 ± 0.04 


* 


Aug 20 08:46:53 


7.920 x 10 3 


PROMPT-5 


Rc 


660.0 


17.52 ± 0.09 


3 


Aug 20 09:03:19 


8.906 x 10 3 


P60 


Rc 


360.0 


17.69 ± 0.04 


* 


Aug 20 09:28:50 


1.044 x 10 4 


P60 


Rc 


360.0 


17.85 ± 0.04 




Aug 20 10:00:44 


1.235 x 10 4 


P60 


Rc 


360.0 


17.78 ± 0.04 


* 


Aug 20 10:38:59 


1.465 x 10 4 


P60 


Rc 


360.0 


17.97 ± 0.04 


=i= 


Aug 20 11:05:54 


1.626 x 10 4 


P60 


Rc 


360.0 


18.01 ±0.04 


* 


Aug 20 11:33:14 


1.790 x 10 4 


P60 


Rc 


360.0 


18.03 ±0.04 


=i= 


Aug 20 12:02:31 


1.966 x 10 4 


P60 


Rc 


360.0 


18.16 ±0.05 


* 


Aug 21 00:37:53 


6.498 x 10 4 


RTT150 


Rc 


900.0 


19.36 ± 0.01 


1 


Aug 21 04:16:53 


7.812 x 10 4 


PROMPT-5 


Rc 


5370.0 


19.94 ±0.31 


3 


Aug 22 00:00:05 


1.491 x 10 5 


RTT150 


Rc 


1800.0 


20.26 ± 0.05 


1 


Aug 22 07:17:03 


1.753 x 10 5 


P60 


Rc 


6840.0 


20.51 ±0.11 


=i= 


Aug 22 23:00:41 


2.319 x 10 5 


RTT150 


Rc 


3900.0 


20.90 ± 0.03 


2 


Aug 23 08:06:11 


2.647 x 10 5 


P60 


Rc 


8400.0 


20.89 ±0.10 


* 


Aug 23 22:18:05 


3.158 x 10 5 


RTT150 


Rc 


2700.0 


21.18 ±0.04 


2 


Aug 24 08:46:32 


3.535 x 10 5 


P60 


Rr 


8400.0 


21.22 ± 0.11 


* 


Aug 25 09:33:11 


4.427 x 10 5 


P60 


Rr 


2880.0 


21.34 ± 0.13 


* 


Aug 26 05:05:39 


5.130 x 10 5 


HET 


Rr 


600.0 


21.57 ± 0.08 


* 


Aug 26 08-28-20 


5 252 x 10 5 


P60 




3600 


21.64 ± 0.12 


* 


Aug 27 08:37:13 


6.121 x 10 5 


P60 


Rc 


4800.0 


21.80 ±0.12 


* 


Aug 27 22:49:53 


6.633 x 10 5 


RTT150 


Rc 


1500.0 


22.02 ±0.10 


4 


Sep 26 01:39:38 


3.179 x 10 6 


HST 


F625W 


800.0 


24.59 ± 0.08 


* 






Rc 




24.55 ± 0.08 


=i= 


Aug 20 06:40:21 


328.0 


P60 


k 


60.0 


14.91 ± 0.02 


=1= 


Aug 20 06:45:19 


626.0 


P60 


k 


60.0 


14.42 ± 0.01 


* 


Aug 20 06:50:06 


913.0 


P60 


k 


60.0 


14.78 ± 0.01 


* 


Aug 20 06:56:37 


1.304 x 10 3 


P60 


k 


120.0 


15.24 ± 0.02 


* 


Aug 20 07:07:15 


1.942 x 10 3 


P60 


k 


120.0 


15.74 ± 0.02 


* 


Aug 20 07:18:03 


2.590 x 10 3 


P60 


k 


120.0 


16.07 ± 0.02 


* 
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TAB LE 8 — Continued 



Mean Observation Date 


'bat 


lelescope 


Filter 


Exposure lime 


Magnitude a 


Keierence 


\ jL\JVjJ U 1 ) 










I vegaj 




Aug 20 07:40:11 


3.918 x 10 3 


P60 


k 


360.0 


16.54 ± 0.03 




Aug 20 08:21:15 


6.382 x 10 3 


P60 


k 


360.0 


17.02 ± 0.05 


* 


Aug 20 08:46:26 


7.893 x 10 3 


P60 


k 


360.0 


17.22 ± 0.04 


* 


Aug 20 08:46:53 


7.920 x 10 3 


PROMPT-3 


k 


1560.0 


17.31 ±0.08 


3 


Aug 20 09:11:47 


9.414 x 10 3 


P60 


k 


360.0 


17.23 ± 0.05 


* 


Aug 20 09:37:27 


1.095 x 10 4 


P60 


k 


360.0 


17.34 ± 0.03 


* 


Aug 20 10:21:01 


1.357 x 10 4 


P60 


k 


360.0 


17.48 ± 0.04 


* 


Aug 20 10:47:55 


1.518 x 10 4 


P60 


k 


360.0 


17.68 ± 0.04 


* 


Aug 20 11:14:54 


1.680 x 10 4 


P60 


k 


360.0 


17.71 ±0.05 




Aug 20 11:42:22 


1.845 x 10 4 


P60 


k 


360.0 


17.72 ± 0.04 


* 


Aug 21 04:04:53 


7.740 x 10 4 


PROMPT-3 


k 


5440.0 


18.33 ±0.11 


3 


Aug 22 00:07:53 


1.496 x 10 5 


RTT150 


k 


1800.0 


19.74 ± 0.07 


1 


Aug 22 07:25:05 


1.758 x 10 s 


P60 


k 


6840.0 


19.97 ±0.12 


* 


Aug 22 23:06:41 


2.323 x 10 s 


RTT150 


k 


3900.0 


20.37 ± 0.05 


2 


Aug 23 08:22:42 


2.657 x 10 5 


P60 


k 


7560.0 


20.48 ±0.12 




Aug 23 22:24:05 


3.162 x 10 s 


RTT150 


k 


2700.0 


20.78 ± 0.09 


2 


Aug 24 09:11:48 


3.550 x 10 5 


P60 


Ir 


8400.0 


20 59 + 1 1 


* 


Aug 25 10:44:39 


4.470 x 10 5 


P60 


Ir 


2400.0 


21.02 ± 0.15 


* 


Aug 26 04:49:33 


5.121 x 10 5 


HET 


Ir 


1200.0 


21.25 ± 0.10 


* 


Aug 26 09:27:24 


5.288 x 10 5 


P60 


Ir 


4800.0 


21.17 ± 0.14 




Aug 27 08:34:02 


6.120 x 10 5 


P60 


Ir 


4440.0 


21.30 ± 0.13 




Sep 26 02:01:58 


3.180 x 10 6 


HST 


F115W 
k 


800.0 


24.32 ± 0.09 
24.27 ± 0.09 


* 
* 


Aug 20 06:41:50 


417.0 


P60 


z' 


60.0 


13.93 ±0.11 




Aug 20 06:46:45 


712.0 


P60 


z' 


60.0 


14.27 ±0.14 


* 


Aug 20 06:51:50 


1.017 X 10 3 


P60 


z' 


60.0 


14.62 ± 0.21 


* 


Aug 20 06:59:16 


1.463 X 10 3 


P60 


z' 


120.0 


15.02 ±0.13 


* 


Aug 20 07:09:56 


2.103 x 10 3 


P60 


z' 


120.0 


15.49 ±0.14 


* 


Aug 20 07:20:45 


2.752 X 10 3 


P60 


z' 


120.0 


15.96 ± 0.20 


* 


Sep 26 03:06:14 


3.184 X 10 6 


HST 


F850LP 


1600.0 


24.09 ± 0.09 


* 



REFERENCES. — * = this work; 1 = IBikmaev et alJ I2005|) : 2 = IKhamitov et alJI2005l) : 3 = IMacLeod & Nvsewandeil 120051) : 4 = IAslan et all 120051) . 

a Errors quoted are lc photometric and instrumental errors summed in quadrature. Galactic extinction (E(B — V) = 0.044; Schleael et al. 1998) has been 
incorporated in the reported magnitudes. 



